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THE FIRST PART OF A MEMOIR ON THE DEVELOPMENT OF THE 
DISTURBING FUNCTION IN THE LUNAR AND PLANETARY 
THEORIES. 


[From the Memoirs of the Royal Astronomical Society, vol. XXVIII. (1866), pp. 187—215. 
Read November 10, 1858.] . i 


THE development, as is well known, depends upon that of the reciprocal of 
the distance of the two planets: and Hansens Memoir “Entwickelung der negativen 
und ungeraden Potenzen der Quadratwurzel der Function 7?+ r° — 2rr’ (cos U cos U’ + 
sin Usin U’ cos J),” Abh. der K. Sáchs. Ges. zu Leipzig, t. 11, pp. 286—376 (1854), 
contains a formula which is truly fundamental, viz. the expression of the coefficient 
of the general term 

n 

Tn cos (JU +7 U’) 
of the development of the reciprocal of the distance as expressed in the above- 
mentioned form, where r, 7” are the radius vectors of the inferior and superior planets 
respectively, and U, U’ are the angular distances from the mutual node. In the 


lunar theory, where the higher powers of are neglected, we have in this manner 


a small number of terms each of which is to be separately developed in multiple 
cosines of the mean anomalies. This can be effected as in the Fundamenta Nova, 
and my “Memoir on the Development of the Disturbing Function in the Lunar 
Theory,” R. Ast. Soc. Mem., t. XXVIL, 1859, [213], which is a mere completion of 
Hansen’s process('). In fact if f, f’ are the true anomalies, and ©, ©’ the distances 


11 take the opportunity of mentioning the memoir of Hansen’s which immediately precedes that above 
referred to, viz. “Entwickelung des Products einer Potenz des Radius Vectors mit dem Sinus oder Cosinus 
Cines Vielfachen der wahren Anomalie in Reihen die nach den Sinussen- oder Cosinussen der Vielfachen der 
Wahren der excentrischen oder mittleren Anomalie fortschreiten,” t. 11., pp. 188—281 (1853). 


www.rcin.org.pl 


320 THE FIRST PART OF A MEMOIR ON THE DEVELOPMENT OF THE [214 


of the pericentres from the mutual node, then we have U=/+6, U=f' +0, 
and the general term is 


ee . “pi ` Wa? YI 
wan os (f+ jf +j6 HIO) 


where r, f are given functions of the mean anomaly g, and 2”, f’ are the like functions. 
of the mean anomaly g’. And the development depends upon those of 


pa ig a po 


sin sin 


which (if we consider as well negative as positive values of the index n) are each 
of the form 


and when the developments of these expressions are known, we obtain at once by 
the mere addition and subtraction of the coefficients of the cosines and sines of the 
different multiples of g and g’, the development of 


d belegt vd D 
pnt cos OU LIT +96 KS O ) 


in the tabular form employed in my memoir just referred to. In the planetary 
theory we must unite together the terms containing the different powers of i SO as: 


to form the entire coefficient D (j, 7’) of cos (jU+jU); if then we write 
r=a(1+2), WT (1+2’), 
and develope the coefficient in powers of a, a we have the general term 
x ae cos (jU+3jU”) 


which admits of development in multiple cosines of the mean anomalies, in the same 
manner precisely as the before-mentioned general term 


== ‘ art 
wari 008 (GU +) U’) 

in the lunar theory. It is proper to remark that this method is really identical with 
that commonly made use of in the planetary theory: the only difference is, that by 
Hansen’s fundamental formula, we have the complete expression of the coefficient 
DO. j) developed in powers of sin or of tan ké, instead of (as in the ordinary methods) 
the first two or three terms of this development. 
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The required development of 


a x“ cos (jU JU) 
or, what is the same thing, 


eer cos (if if +j6 HIO) 


depends on the developments of 
cos . jh OOS o is 
on if, eene 


These are functions of the same form, and we may consider only 


cos. 


sind 


The value of z is (= - 1) and we could of course calculate 


r * COs . 
ec ei, 

by the methods of the Fundamenta Nova or the memoir of Hansen’s referred to in 
the foot-note. But if we write f=g+y (y is the equation of the centre), then the 
required expressions depend on 


which are actually calculated as far as e” for a=0, 1, 2...7, and j an undetermined 
symbol, by Leverrier in the Annales de lObservatoire de Paris, t. I. pp. 346—348 
(1855). Hence, by the mere substitution, in Leverrier’s formula, of the numerical 
values of j and f, and by the addition and subtraction of the coefficients of the 
cosines or sines of the different multiples of g and g’, we may obtain the develop- 
ment of 


aea cos (jf + jf’ HIO HIT) 
in a tabular form similar to that employed in my memoir already referred to. 


I have thought it desirable to put together the various results above referred to, 
and to investigate by a different process the expression for Hansen’s coefficient, which 
in his memoir is obtained by means of a long series of transformations which it is 
not very easy to follow, and is not exhibited in quite the most simple form. And 
this is what I have done in the present first part of a memoir on the development 
of the disturbing function. My object has been to exhibit, in as complete a form as 
possible, the preliminary development in multiple -cosines of the true anomalies; and 
to indicate the process of the ulterior development in multiple cosines of the mean 
anomalies, 
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E 
Let M be the inferior, W’ the superior of the two planets (in the lunar theory 


M is the moon, YY the sun) and let the quantities relating to the two bodies be, 
for M, ; 


r, the radius-vector, 

U, the distance from node, 

©, the distance of pericentre from node, 
f, the true anomaly, 

g, the mean anomaly, 

a, the mean distance, 

e, the eccentricity, 


and for M' the accented letters, 7’, Ze, in the like significations. 
The node referred to is the ascending node of the orbit of M upon that of W, 
and I write also 
®, the inclination of the orbit of M to that of W, 
y, =sin¿0. 
The disturbing function is in the first instance given as a function of r, 7’, H, 
where 
cos H = cos U cos U’ + sin U sin U’ cos O, 
that is, as a function of r, 7’, U, U”, ®, or, what is the same thing, of r, 7’, U, U’, a 
And the preliminary development is a development in multiple cosines of U, U’. We 
have then 


U =f +0, 

U=f' +o’, 
and finally e , 
=a elqr(e, y), 


f = elta(e. g), 
r” =a’elar(e’, y”), 
F = elta(e, gl 


and the ulterior development is a development in multiple cosines of g, y”, D Di, the 
coefficients involving, as before, 7, and also a, e, a’, e. But as usual it is not 
attempted to carry the development further, by introducing in the coefficients in place 
of the relative quantities UG, O’, 9, the remaining elements of the two orbits, which, 
if it were necessary to use them, would be 


for W, 
0, the longitude of node, 
o, the departure of node, 
$, the inclination, 
w, the departure of pericentre, 
and for W, the like accented quantities, the orbit of reference being any fixed or 
moveable orbit whatever. 
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IL 


The expression for the disturbing function on M- (that is, when the superior 
planet disturbs the inferior) is 


ml 1 
r Vr? + r?2—2rr' cos H 


and that for the disturbing function on W (that is, when the inferior planet disturbs 


the superior) is 
r cos H 1 
M + — La 
| yA Nri + 7? — rr =H 


where the disturbing function is taken with Lagrange’s sign (=—R, if R be the 
disturbing function of the Mécanique Céles te). 


But we may in the first instance consider the development of the reciprocal of 
the distance of the two planets 
KEE 1 
Nr? + 12 — 2rr’ cos H ` 


The preceding expression for cos H may be written 
cos H = cos U cos U’ +sin Usin U’ (1 — 277), 


or in either of the two forms 
cos H = cos (U — U’) — 27? sin U sin V’, 
cos H = (1 — y?) cos (U — U’) + n? cos (U + U’). 


Now imagine the function developed in ascending powers of 5, the coefficient of 


qe e S ; S 
den will contain cos” H, cos” H,... to cos H or 1 according as n is even or odd; and 


if we then substitute for cos H the last given expression, and express the different 
powers of cosH in multiple cosines of U—U” and U+U’, and make the final 
expression contain the cosines of opposite arguments each with the same coefficient, 
it is easy to see that the form of the general term is 


an On (J, J) cos(jU+5U, 
where j, 7’, each of them extend through the values n, n—2,...—n, and where 


Cn(—j, —J) = Cn (j, J). 


; 41—2 
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Thus in particular for n=0, 1, 2, 3, the combinations (j, Jj) belonging to the several 


arguments are, 


But as the coefficients C satisfy the condition C(—j, —j’)=C(j, J) the two terms 
C(j, J)cos(jU JU) and C(=j, —j)cos(-jU—j'U”) are equal to each other, and they 
may be combined together into a single term. The general term may consequently be 
` written 
an 20a (j, 5) cos (jU+5U, 

where j has only the values n, n—2,...1 or 0, and Jj has as before the values 
n, n—2,...—n; except (which occurs only when n is even) for 7=0, when j’ has only 
the values n, n—2,...0: and in the particular case 7=j’=0, the last-mentioned 
expression for the general term must be multiplied by 4, or, what is the same thing, 
the factor 2 must be omitted. In particular for n=0, 1, 2,.3, 4, the combinations 
(j, J) belonging to the several arguments are 


SSS LO 


EPT A eta lä, ZO, ¿2 a aes ee län 4/2, 4|0, 4 
biet 2, OJO, O EC ck es WA 
2,—2 ade E it ed 4, 0/2, O e TA 
P 30-3 | Zë? 4,-2|2,-2 Kc 
4,—4)| 2,3—4 


I remark that in a series K (j, 7’) yor jU JU, where each argument occurs positively 
and negatively, and K(—j, —j’)=+K(j, Jj) according as the series is one of cosines 
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or of sines, we may say that the discrete general term is K(j, Jj) yA (GU+3U”, or 

that K(j, Jj) is the discrete coefficient of the cosine or sine, but if we unite 

together the terms with opposite arguments so as to form the general term 


Ea, $) = (jU + 7'U’), then this may be called the concrete general term, and 


2K (j, J) may be said to be the concrete coefficient of the cosine or sine. In a sine 
series the term corresponding to the argument zero vanishes; in a cosine series this 
is not in general the case, and the concrete term corresponding to the argument zero 
must be multiplied by A 


Returning to the question in hand, from the symmetry of the expression to be 
developed, we have 


Cn (j; JI =Cy, Gf Y); 


and it follows that the only coefficients which need be calculated are those for which 
j is not negative, and not less in absolute magnitude than 7’; the remaining coefficients 
are respectively equal to coefficients which satisfy these conditions. Thus for n=0, 1, 2, 3, 4, 
the combinations (j, 7’) corresponding to the coefficients in question are, 


0, “Gis 14g 


and we have, for instance, O (2, 4)=C(4, 2), C(- 2, —4)=C (2, 4)=C (4, 2), &. Under 
the preceding restriction, viz. j not negative, and not less in absolute magnitude 
than J, the expression for C,(j, 7’) (deduced from the formule of Hansen’s Memoir) 
is as follows; viz. putting as usual Ilz=1.2.3...7, and also Il, (u—4)=34.3.3...(x—J), 
and representing the hypergeometric series 

a.8 a.a+1.8.8+1 


1 + a+ 


A O a g P 
JE? at E 


by F(a, B, y, 2), we have 


 — 2G (n 49) - PG nn 
CO N= Hy @—p Ud —j) G+) 


x nit (1 PAD F(—n+j, n+ 941, ¡+ +1, 7) 
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which, it is to be noticed, is a rational and integral function of y, the highest power 
being y”, and the lowest power or order of the coefficient being ai I reserve the 
demonstration of this formula for a separate paragraph. 


Let the discrete general term involving cos(jU + ’U’) be represented by 
D (j, J) cos (GU zU 
we have, in like manner as for the coefficients C, D(-j, —7) =D (j, j), D(F, ss DO, J), 
and consequently D(j, J) will be known, if we know its value when the before- 
mentioned conditions are satisfied; viz, if 7 be not negative and not less in absolute 


magnitude than 7’; and collecting together the different terms which involve the cosine 
in question, we find at once, the conditions being satisfied, 


E at io De"? 
D(j, I) = T C3(9, J) + pits Ciro (J, J) + &e.... 


so that the value of DO. 7’) is known. A transformation of this expression will be 
given in the sequel. 


TIT. 
The concrete general term is 
2D (j, $) cos (jU +j U’); 
in particular the concrete terms involving the arguments U—U’, U + U’, are 
2D (1, — 1) cos (U — U’), 
2D (1, 1)cos(U + U”, 


or, substituting for the D coefficients their values, the terms are 


2 1 C, (1, - 1+7,0,(1,-1)+ sa) cos (U — U’), 


son, geben, D+} eos(U+ 0) 


So far I have considered the reciprocal of the radius vector; but if we consider, 
instead, the disturbing functions, the only difference will be that we must add the term 


i 
— ¿2005 H, 
or 


Pi 
— za COS H, 
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according as the superior planet disturbs the inferior one, or the inferior planet the 
superior one. The value of cos H is 


= (1 —7*)cos(U—U’)+ op cos (U + U’), 
= 20,(1,—1)cos(U— U’)+20,(1, 1) cos(U + U’), 
observing that 


C1-D=3+(01-7) Gd, 1)=)7; 


and the terms to be added are consequently, when the superior planet disturbs the 
inferior, 


2 C,(1,— 1) cos (U — U’) 


2 Ze C,(1, 1)cos(U+0", 


the effect of which is simply to destroy the same terms contained with the apo 
sign in the reciprocal of the distance; 


and when the inferior planet disturbs the superior one, the terms to be added are 
— GO, — 1) cos (U — U’) 
— 2 30 (1, 1)cos(U + U”, 


which are not equal to any terms in the reciprocal ot "he distance. We may write 
as follows : 
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Reciprocal of Distance is, 


Disturbing function + Mass of Disturbing Planet is, 


U +U 


D (o, 0) 
+2 D (1, —1) 


O| -230 -1) 


EST te LA EA 


+ 

D 

o 

2 

N 

w 
w ® o NN 
LENNON 


+ 
N 
SÉ 
9 
m 
== 

bd QQ Mm QQ 

[lp 

ka WF H 


+ 
N 
S 
a 
E 
o 
ks 
=D &. bk b Doe e 
| 
Sab Rb O RE Em em 


Ps Only to be inserted for the disturbing function when the superior planet disturbs the inferior, and 
having the effect of destroying portions of the coefficients of the next preceding terms. 


2 Only to be inserted for the disturbing function when the inferior planet disturbs the superior one. 
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For the lunar theory, to the extent to which it is necessary to carry the develop- 
ment, and to which it is carried in Hansen’s Fundamenta Nova, and my memoir 
before referred to, we might simply write, 


Disturbing function + Sun’s Mass is 


(4) 
(5) 
(7) 
(9) 
(10) 


I 
e C, (o, 


g? 
= 73 Co (o, 


qa 
+27 Cs, — 


e de 


qa 
13 Ca (2, 


+23 
+27 bala 
+250,(0, 
+25 0,(2, 
+2 Set 
+2 5,053, 
e bat 


qu 
+ 2 7 Os (3, 


Ae afin 


U +U 
as ay 
o O 
fe) O 
A =P E 
I I 
71 AA 
2 (e) 
o 2 
2 2 
lala 
3 -1I 
I -3 
3 I 
CR 8 


where the prefixed numbers are those of Hansen e ten parts Q,, Q... Quo; or omitting 
the first term which depends only on the sun’s radius vector, and the last two terms 
which are ultimately neglected, and substituting for the coefficients C their values we 


have, 
G. TIE 
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Disturbing function + Sun's Mass 1s 


(1) 

(2) 

(3) 

(4) $y- dr S o 2 
(5) SI y 2 2 
(| ie nget [ae 
(7) 5 Bei Bel 323 
(8) tee za SE $ 
(9) ere 2 gts 
(10) dia D cre 3 


If the same form be adopted for the planetary theory, the expressions for the 
leading coefficients will be, 


C, (0, 0), 


2p 


D(0,0)= 5 1 


t- Gert, Gei 

yé 
+ ao 90 — 1407+ 707!) 
He (1 — 42 n? + 420 nt — 1680 ag + 3150 n°...) 
+ m 4225 (1 — 72 n? + 1260 nt — 9240 n° — 34650 n°...) 


YA MOPD AL 


y/2p+1 Un Un 
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a, 1), DA. (=. T 


ya 


+h An (-107 + 157) 


331 


n°); 


n°); 


5 +7, Be =n) (1-28 9? +168 nt — 336 9° + 210 7°) 
7 
7 += ALS, (1 — 9?) (1 — 547 + 675 an — 8100 ap + 7425 7°) 
rev) TL, (p +4) Th (p= A 2). 
2p +1 + erg TO o ao aa een 2p +3, 1, 7°); 
C. (1, 1), DU, us A Ae 
éi H zb 0 m2 4 
3 La ¿TU-YN+ Bei 
5 +5, $38 1° (1 — 38 n° + 289-188) 
7 + Ta BAP of (1 — 18 7° +230 mf 33017) 
reno ¿TL (p+4) Ih (p+3) ; 
2p+1 + em 2 Tp Tp Y F(— 2p, 2p +3, 3, n°); 
r? i 
C: (2, — 2), D (2, —2)= rs 3 (1 — 2? 
Ku 
H + $ (1 — 17°)? (1-14 + 28 9!) 
6 +5 aie (1 — n°)? (136 9? + 270 9! —660 + 495 1) 
; eg Th (p +), (8-8 br, éi 
2p tan TE pI) pal) 7 FC P+ % 2 +3, 1, 
C e 2 2 
E 2 2 14 m2 4 
4 ta PA-P) A- n 4) 
6 + $88 at (1 — af) (1 = 12+ 45 94 — 66 9... 
y TL, (p+) th q 4) Y e FH- © 3 3 


42—2 
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Cr (0, 2)=Cn (2, 9), D(0, 2)=D(2, 0); 
C. EE 
2 (2, 2), D (2, 2)= pe $7 
4 r 5 pi (1 — 14 m + 28 nt 
ta TL E + T) 
6 ” 367 4(] — 326 n2? + 18 nt 
Le BP n (1 — 3E n? + 1891) 
re 1, (p +4) Th (p +d) nc, 
2p + ven Hieft 7" 2p +2, 2p+3, 5, n°); 
r 
C: (3, — 3), 28-12 Ce: elegy 
5 ” a 1 —2P(1 2 4 
+ 77 tes (1 — n°) (1 — 180? + 457°) 


rr Th, (p+ $) Uh (p — §) 


par TI (p +2) I (p—1) (LP (= 2p +2, 2p + 5, Lais 
r? 

C; (3, — 1), D(3, —1)= y 4 15 7? (1- 7?) 

5 +5, Mal — PY (1 pat 

ræt | TI, (p +4) Ih (p-4) s 

2p +1 + aM EC wett tb än +2, 2p+5, 3, 7°); 

DH p 

CO, (1, —3)=C,(8, — 1) DE, ss DS 

C ii 4 2 

5 +5, BEMA- E Aere 

; r#+ TL (p +3) IL, (p+ ; Ya 

ST a E F(-2p+2, 2+5, 5,7); 

Cali, Biz Dat, Ir Dil, D= DB; "D: 

O, (3, 3) Dog tae 

3 NG; ` D = va Tg” 

5 za 5 nê (1 418 p? 

SET dE a...) 

2p-+1 "ef 

PEI E ty PED WAA ar F(— 2p +2, 2p+5, 7, 1); Be. 


pra de 1 (p=1) Mp1) 
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IV. 


But for the planetary theory it is more suitable to arrange the expression of 
D(j, J') according to the powers of a We have, under the before-mentioned conditions, 
j not negative and not less in absolute magnitude than 7, 


LF sf iter haji 
DX}; J) = > am Cra (3 J) 


ad pa Mes 
= 3 (5) Cira (J, J’) 
where A extends from 0 to œ; and, writing in the expression of C,(j, 7’), j +2) 
for n, we find 


or, substituting for the hypergeometric series its value, 
tie TIA IES (j=) (+7) (Oy [5+3+ OY 
e ST EE ET) 42004 0) ae: 
e DAI- G47 Oat UCA- AUA 


from 9-0 to O=w. Substituting in the numerator for 1121, 2A TIA T, (à -— +4), and 
in the denominator for I (2j +2), 2*A II (+A), (¿+A-—%3), and reducing, this 
becomes 


o ` ny (=) n? 1 TL, G(G+ +2 DI (9 + 22, + 0) TL (A — 4) 
nalh D)= Gm WF MEG NINIGINMECAO > 


and substituting this expression in D(7, 7’), we find 


+ HUWA CE 
DG, j) = gy PI A-H Ze ong NT 


a (© GJ rh H (j++ 0) TL H CT 


IL (4 (j 7) +d) 1 (¡+ 120-8) 
from A=0 to A=% and 0-0 to 9-0; which is the required expression. It is to 
be remarked that the series in (2) which multiplies 9? is not in general expressible 


as a hypergeometric series. If, however, we attend only to the leading term in y, 
or write 0= 0, we find for D (j, 7’) the value 


ee ga bing 79 1 
mia AA KU ar ai 
23 AEN Geen riit T1(9 +7/) 


a (7 (4 (fj +7) +%— Y) (A+A (A— 3) CT) 
7 


IEG (j—J) +») (G+) 2d 
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which may be simplified by putting in the numerator for II (2j + 2), 
Zit TI (j +A) Th (j +A — $), 
and in the denominator for II (j +j’), 
PINK +7) Th G G+H) — 9), 
and for 1121, 2A III (à— 4) We thus obtain the value 
1 
T, “G +) UG G+) —D 


EE 
MEG) +4) Un (7) > 
which is equal to 


mp 
nit (1— pain Ze a EGER 


E 
(39-7) + APP EE 
or finally we have, as regards the leading term in 7, 


i IL, (j — 3) 
D(j, said Ou? d 


FIGO HI hita AA, 5). 
I remark that if in general 
1272 (72 + r? — 2rr' cos Y) *+=23R,¿cos1Y, Ki Re 
then, writing $(j—7’) for i and 4(7+,’) for z, we have 
Re = Fi EG DA pT EGA th ith KA Ra YA $) 
and the last-mentioned expression for D(j, Jj), as regards the leading term in 7, 
becomes therefore 


SN 34 yu MEGAN) pus 
DG Den TADB A 


and we have in general 


DG, =n (1 — pps Ke: ib Gei m > De go iris d . 


As a verification, it may be noticed that for 7+ j’=0 we have DU. — j) = (1 — nY 
(Rj + terms in y), and for »=0, D(j, — j)=Rj4, which is right, since the two sides 
each denote the coefficient of cosjY in (r?+r2-—2rr'cos 9) There is reason to 
believe that the expression for D (j, j’) might be further reduced so as to obtain in 
a convenient form the coefficients of the successive powers of y; but I have not yet 
accomplished this. 
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Y 


Considering now D(j, f) as a given function of r and 7’, we must in the 
planetary theory write r=a(1+~2), r= (1+2), and develope in powers of z and 


7. The general term is, of course, 


Mama dede DO Die 


where DC. j’) is what D (j, Jj) becomes when a, a’ are substituted for r, 7’; and, 


writing f+, CL for U, U’, we see that in the planetary theory the terms to 
be developed are of the form 


xx” cos (Jf +7 f' +jG6+ 70’), 
while, from what has preceded, the form for the lunar theory is 


rar” cos (jf +7f' +j64+70), 
(x is always =—n—1, except in the terms = cos (U +U d 


the values which have to be substituted being 


r =a elgr(e, y), æ =elqr (e, g)—1, 
f= elta (eg); 
r =a elqr (e, y”), x = elqr (e, y) —1, 
f'=  elta (e, y”). 
Suppose that in the former case the developments of 
arcos jf, arsinjf, event, es sin jf’ 


and in the latter case the developments of 


$ 


s\n n WEN nv 
(2) cos jf, (2) sin jf, (E) Ge G sin ZC" 


a 
are 

> [cos]'cosig, X [sin] sin ig, > leo oos do, 2 [sin] sin do, 
where the summations extend from i or “=—o to o, and where the coefficients 


[cos], [sin]? satisfy 
[cos] += [cos], [sin] +=- [sin], 


and in like manner the coefficients [cos]’, [sin]’ satisfy 
[cos]-* =[cos]",  [sin]-* =—[sin]’. 


It is to be observed that (cos), [sin]? are functions of e, and [cos]*, [sin]* 
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functions of e. the accents to the indices 7 and Y are sufficient to indicate this. 
Hence observing that 


> loost cos ig . 2 [cos] cos y = 22 [cos]'[cos]" cos (ig +79’), 
3 [cos] cos ig . X [sin] sing = EE [cos]! [sin] sin (ig + 7g’), 
> [sin]' sin ig . = [cos]" cosz’g’= 22 [sin]! [cos] sin (ig 4797), 
X [sin]! sin ig . 2 sin)” sin Co = — 22 [sin]! [sin]’ cos (tg +19”), 


S mM fy! Nn 
we have for the products of oa, or as the case may be KR (=) into 


cos (jf+jf’), the values 22 (ege [cos]* + [sin]' [sin]*) cos (19 + 19”), 
sin (yf LIZ) > 22, (amb [cos]” + [cos] [sin]") sin (19 + t'g’), 
and thence observing that 
cos (jO +70’) ÈX Pi" cos (ig + ig) = 22 DT cos (ig + Uy" +30 +30"), 
— sin (¡UT +0") EX Qb” sin (ig +19”) = TE Qi” cos (ig IO +30”), 
provided only that P-5 =P", Q- =- Qb", we find for the product of Ger op 
as the case may be CT ER into cos(jf+jf' +jO+ jG’) the expression 


=> ([cos]! (cos)? + [sin] [sin]* + [sin] [cos]” + (cos) [sin]”) cos (ig + t'g’ +76 +30") 
or, finally, the expression 
23 ([cos]' + [sin]*) ([cos]* + [sin]”) cos (ig + vg’ +76 +30") 


which is the required development of the general term in multiple cosines of the 
mean anomalies. 


Vib, 


Investigation of the coefficient Cn (J, 7’). 


It is possible that there might be some advantage in developing in the first 
instance according to the powers of cos H, a process which, as it has been seen, 
leads very readily to the form of the general term; but the mode which I have 
adopted is to develope in the first instance according to the powers of y. I put 
therefore 

cos H =cos(U — U’) — 27? sin U sin U’, 


and we have then for the reciprocal of the distance, 
{7? +1 — 2rr cos (U — U’) —rr' (- 4 sin U sin U’) y) 
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which is to be developed in ascending powers of z and we have for the discrete 


general term of the development 
a H Hi a Va 7 
nt Cr (j, J’) cos (JU +7 U’); 


which is the definition of OO. 7’), the coefficient the value of which is sought for. 
It has been seen that j, Jj, are each of them of the same parity with n (even or 
odd according as n is even or odd), and it has been seen also that it is sufficient 
to consider the case where j is not negative and not inferior in absolute magnitude 
to 7’. 


Expanding in powers of y, and putting as before 
Mx=1.2.3...%, TI, (1—-4)=3.3.8 ... (@— 4) 


the general term is 
Lin pay? ((1 + 72 — 2rr cos (U — U”) 793 (- 4 sin U sin UY y” 
where œ extends from 0 to o. 


The factor (—4sin U sin U’) consists of a series of multiple cosines, and as usual 
it is assumed that the cosines to opposite arguments are made to occur with equal 
coefficients. The form of the general term is cos(AU +w U’), where A, A have each 
of them the values o 2-2, w—4,...—w, that is, A, A are each of them of the 
same parity with æ. Hence j, 7’ being as before of the same parity with each other, 
and Y being even or odd according as j, j and æ are of the same parity or of 
opposite parities, the development of (—4sin U sin Ur will contain a series of terms 
cos [(j +9) U+(3' — Y) U], which (since the other factor contains only multiple cosines 
of U-— U’) are the only terms which give rise to a term cos(jU +7 U’). I represent 
the discrete term of (— 4 sin U sin U”* which contains the before-mentioned argument by 


MS cos [(¿ +3) U+(3- 9) U]. 


On the before-mentioned assumption, 7 not negative and not less in absolute 
magnitude than J, we have j +f and ¿—7, each of them not negative. We must 
have j+% $a and j —S}a, that is ¥$(e—j) and + — (x -— 7’), consequently 
?2—3j4—(r—J) or 2at¢j+7’. And this relation being assumed, Y extends from the 
inferior limit —(«—j’) to the superior limit (e —7) by steps of two units, the extreme 
terms being . 

sie Ma? cos[(j +j — 2) Defi 


$= («-j), Mè cos [2U+(j+3' —2) VU) 


the coefficient of U increasing from j +f —e to æ, and that of U’ diminishing from 
æ to 9+ 7’—a« by steps of two units. 
Cs EE. 43 
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. . . y . 
Now expanding in ascending powers of za write 


rey! ((72 +r —2rr' cos (U — U’)}-*4 = ZR cosi (U — U’) 
where 1 extends from —o to +o and R= R+; so that R,'cosi(U—U’) is the 
discrete general term of the development. The term R,*cos$(U—U") in combina- 
tion with the term M,>cos[(j7+%) U+(j’—9) U’] gives rise to M,*R,* cos (JU JU), 
and restoring the multiplier which has been omitted, and giving to Y and æ the 
different admissible values, we find for the discrete general term containing cos(j7U + 'U’) 
the value 


z EE-D) ya SO Re?) cos (JU + JU”) 


where 3 extends from the inferior limit Y =—(#—,’) to the superior limit $=#—J, by 
steps of two units, and æ extends from #=4(j+ 7’) to e=«. The portion of this 


mn 
containing en cos (JU + 'U’) is 


n+l 
sani On (js J) cos (JU JU) 
and we have therefore 


Cr (j, J) = coeff. a SS Kee Up 2 ONS RAA 


eh at 
There is some speciality in the case 7+ j’=0, but the result just obtained subsists 
without variation. To find M,* we have 


M,* = Discrete coeff. cos [(j +9) U+(3'— 9) U] in (— 4sin U sin U’)*, 
and putting sin U= e (v — >) siny =1 (v - >) where i= —1 we have 
2 uy” 2% v ; 


(— 4 sin U sin UY? = (v- T (v — =) 


and the function on the right hand contains the term 
Ila Hz 
E ai E HH viti Tea, 
OM TFT @=f) ege 
or putting 
2—2f=j7+3, 2— 2231 —9, 
IEEE iC Sh F =H(0—J +9, 
which give integer values for f, f’, since Y is even or odd according as j, j, and 


7 are of the same parity or of opposite parities, and replacing vi v’-¥ by the half 
of its value 2 cos[(j+%)U+(j’—9) U’], the term is 


and therefore 


(HIH EAE a A ENSEM re R.. ESTR 
113 (2-3 IF (@ +7 +9) Tb (@ -j +8) 114 (@ +79) 


x cos [(j + Y) U+(3' — 5) U], 
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and consequently, 


Ie Ie e 
13 (2-9) 13 (e +j +3) UY (2-7 +9) 03 (2479) 


which is the expression for M,?. 


Mè =(-)? 39) 


The expression for R, is found in a similar manner, viz, by substituting for 
cos {U — U’) its exponential expression, by which means the function 


(1? + 7°? — 2rr’ cos (U — TA 


breaks up into a pair of factors, each of which can be separately expanded; the result, 
1 being positive, or zero, is 


RS e II, (1+14+m-3) TL, Ge im 3) Sr 
WË TL (i+ m) Il, (@ — 4) Um TL, (@— $) E 


from m=0 to m=oo: this may also be written 


petit TT, (+11) ; : ya 

os a 
die = gatita DU, (a— 4) F(e+i+}, sti i+1,5); 
writing now D for 2, and +9+2m=n, that is m<3(n—z—Y) (n is of the same 
parity with j, Jj, and % is” even or odd according as j, 7’ and e are of the same 
parity or of opposite parities, m is therefore, as it should be, an integer), R,* contains 
the term 


+. Matt) a) iG @Te—s)~ 9) _ 
en TI} (n— +9) IL, (2-4) I$ (n—a@—-9) IL, (2 Ai? 
en Ke 


if for shortness 


gs= Ns Gmte+S)-)IhGM+e—3)—9) 
* Tit (n—@ +9) Th @— 3) th} (n—2<—9) UL @—-$)’ 

a formula which I assume to subsist as well for negative as for positive values of 9, 
so that Ky *=XK,?. It is to be noticed that if x>nm, then either n—x+3 or n—x-9% 
vanishes, and we have therefore K,*= 


Substituting for R, its value we have 


Cali, I= s HL EH 3) 2 > LUS KA 
where MS KA have the values already obtained, and as before Y extends to 
Y=-— (%— Jj) to Y=(%-—5) by steps of two units, and æ (since K,* vanishes for «<n) 
extends from #=4(j+j') to =m. 
43—2 
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To simplify; write e=3(j+j)+u; S=-(2—J)+2s, =(@—7’) — 2t, so that s+t=u, 
s and ¢ being any integer values (zero not excluded) which satisfy this relation, and 
lastly u extends from u=0 to u=n—}(j +j). We have 


corro KULI ua) 
"frëen "P 


Häss (y UO (j +j) tu) IG (J +7’) +u) 
ý MEI (3 (j +7’) +8) Usj +7) +0)’ 


TG (+ +8 — 4) Il, (3 (0+))+4- 3) 
~ 1G @=j)-H1L SG +7) + 4-3) WE @-J)—-8) ThE G+) +03)” 


and substituting these values, and Bag that the result contains a factor 
1 Qj+j'+u 
JT EIA A EE A oN hi h b e l ced b EE th e 
TEG+H) OO ee TTT GT BH) 
result is 
(-—)” 92u nr 
“IG +i + Za) 


se CED 1 Aa SACH 
MGHS MGH +t Ut ` 1GM-)-s) UG (n—J)—t) 


On (j, J) = 90% ai 2, 


which is easily transformed into 


an gar Au LEGO LEOD- 
OnG, N=FF AA UG © 


if for shortness S = 


s OP BANANA mf FEGH +u +14 n-i) 
Y së" [s] [t] 


s+t=u, u=0 to u=n—}4(j +9); 
the value of the sum S is 
Suma Ft nt jt, j+j +1, ei 
and we have consequently 


fin - Gebeier Ser 


x ali (L—y (an n+j+l, j+] +1, 1); 


GIE Ee 


the required expression for C, (3, j). It only remains to prove the formula for A 
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For this purpose, observing the equation s+¢= wu, I form the equations 


92u e? 1 
[+e + Quy GANADA +31 


[EG +s) +e}? 


[AG +j) + u= E +7) +4)" 


' ECF HG) + ah? 
and thence 


ES Ee = (+j) + uP BG ETGEN 


-IBOR Eu 1 
BFI) Eum WUEMA 


and we then have (putting also (—)” y" = (—)* 78+), 


aa GOTTA FBC- TEEN 


“BG+p)+u 31" TEG T 
(SG +7) +e) ; 


and it is proper to remark that the summation as regards u may be continued 
indefinitely; for, if w=s+t be >n—4(7+ 7’), then one at least of the relations 
s>4(n—-j), t>4(n—j) must hold good, and at least one of the factorials [4 (n —3)F, 
[$ (n—¿)! will vanish, The two factors in the second sum are the general terms of 
two hypergeometric series. In fact, if we put 


$(n+j)+4 =£, 
(ntj) =a, 
E(J HI) +=, 
t(ntj)+e=€-@, 
$(-ntj) =e-B, 


and if, for shortness, we use (aj* to denote the factorial a(a+1)...4a+s—1) of the 
increment positive unity, then we have 


_y fetay weight le ale Bl. 
orn tae SA feta” (Der Er " 


the two hypergeometric series being consequently 


and therefore 


| aoe Br). Cee BOE |! T 

F(a, Bi e+4, a) Gë BEE ai KA EE + &e., 
e e—a.e—B e—a.e—at+tl.e—f.e—B+1 e 
F(e—a, e— B, e+4, A E SS DN L. > nn ni + &e, 
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I assume the truth of the following remarkable theorem, [see 211] viz.: 
“The series formed from the product of the two hypergeometric series, 
F (a, B, e+4, n°), F(e—a, e—B, e+4, 7°), 


e+} e+4.€+3 


by multiplying. the successive terms of the product by 1, ena oe Are, respec- 


tively, is 
= (1 -— petro F(2a, 28, Ze 1)” 
Hence, observing that the general term of S is formed precisely in the manner in 
question, we have 
S= (1— të F (2a, 28, 26, 7) 
or, substituting for a, 8, e their values 
S=(1-PY4 Fi ntj nt+j+l,jt+7 +, 7) 
which is the required value. 
It is clear that a, 8 are interchangeable with e—a, e— 8; that is, we have 
S =(1-— nete: F (La, 28, 2e, y) 
= (1) "Pf F(2e—2a, 2e—28, 2e, y) 
or 
F (2e — 2a, 2e— 29, 2e, y) =(1 — yy +=" F (2a, 28, 2e, n°), 
which is a known property of hypergeometric series. The form 
SIP, ai, jaj, 1) 
is obviously less convenient than the one above mentioned, since the new expression 
is encumbered by a denominator (1 — 407, which really divides out, the finite hyper- 
geometric series containing as a factor the square of such denominator. I have only 


noticed this for the verification which it affords by showing that j, j may be inter- 
changed. 


KE 


The above expression of C,(j, J) is not, in its actual form, given in Hansen s 
memoir. The comparison with Hansens formule is as follows :—The formula (36), 
p. 329, is 

C (n—2f, —(n — 2f — 29)) 

= H cos” 4 J tanvi J F(-(2n-2f-29) — 2, 29 +1, — tan? 4 J) 
where 


Sie II (2n — 2f ) II (2f + 29) 
FT (nf) FM (n—f—g) M(F+g) 129” 
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and the formula (37), p. 330, is 
F (- (2n—2f — 29), — 2f, 2g +1, —tan?$ J) 
= cos 9+4+49 4 J F(-(2n-—2f- 29), 2f+ 29 +1, 29 +1, sim? 3J). 


Combining these, and putting sin4J =n, we have 


C (n— 2f, — (n— 2f— 29)) 
= Hipo (1 — py reto F(—(2n—2f—29), +29 +1, 29 +1, oh 
and then, putting 
! Wich, alte ei 
—n+2f+ 29 =), 


and for OO, j) =0 (3, 7’), writing Ca (j, 7’), we have 


Ca (js F)= Hy (L—p POD F(a tj, nj EL, iti +, ai 
where 
I (n +7 11 (n +3) 


H= anjara- Niin a- iG’ 


or, finally, since 
O (n +j) = 2H I} (n +9) Th G +j) - 9), 
N, (n +j) = 27 03 (n +7’) 11 (1+))- 4), 
we find 


y Z9 Th (+) DM, (+9) ~ 9) 
OO D=TE GP UEO NAGE 


x pH (1D F(—n4+j, n+j+l, + +1, n”, 


which is the formula of the present memoir. 
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